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Abstract 

The Dicke model extended to two bosons of different frequencies or equivalent generalized Jahn- 
Teller lattice model are shown to exhibit a spontaneous quantum phase transition between the 
polaron- modified "quasi-normal" and squeezed "radiation" phase with the transition point depen- 
dent on the frequencies. In a finite lattice a mixed domain of coexistence of the quasi-normal and 
modified radiation phase is created within the quasi-normal phase domain. There occurs a field- 
directed oscillation-assisted tunneling (hopping) . The field is driven by simultaneous squeezing and 
polaron-dressing of the collective boson level mode due to the additional boson mode. In a finite 
lattice in the radiation domain there occurs a sequence of local tunnelings (oscillations) between 
two minima of a local potential weakly coupled to two assisting oscillations. The "radiation" phase 
reveals itself as an almost ideal instanton-anti-instanton gas phase. The correlations among the 
energy levels mediated by the additional mode in the mixed domain considerably reduce the level 
repulsions. As a consequence, the Wigner level spacing probability distribution of the two-boson 
Dicke model is non-universally reduced from the Wigner to the semi-Poisson and asymptotically to 
the Poisson distribution of level spacings. The correlations cause a suppression of the coherence of 
the radiation phase as finite-size effect. Possible applications of the present theory are suggested. 
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I. INTRODUCTION 



The class of nonintegrable pseudospin-boson lattice models with one 



flflflflaaa 



ll| gained a long-term interest 



or two boson modes of different local symmetry 

Bay 

in condensed matter physics and quantum optics. They exhibit a rich variety of interesting 
statistical properties, e.g., quantum and thermodynamic phase transitions and quantum 
chaos in energy spectra revealed in spectral statistics. 

In the Dicke model [5| the quantum spontaneous phase transition was found from the 
effectively unexcited "normal" phase to the "super-radiant" phase, a macroscopically excited 
and highly collective state jf], [7]. A collective spontaneous emission of coherent radiation 
is due to the cooperative interaction of a large number of two-level atoms all excited to 
the upper state (Dicke maser). The transition from quasi-integrability to quantum chaos 
is given by a localization-delocalization transition as a precursor of the phase transition 
in which the ground state wave function bifurcates into a macroscopic superposition for 
any finite number of atoms. Theoretical and experimental research of the superradiance 
or superfluorescence of photons and sp ontaneous emission of phonons in various variants 
" tl 



re Dicke-like collective system 
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12, 



14], ll5|, [16|) still remains in a focus of interest 



If a dephasing mechanism suppresses the strict coherence, processes based on collective 
spontaneous emission will reveal incoherent nature depending on the strength of the dephas- 
ing and some peculiar effects result: Luminescence spectra of rare gas solids, alkali halides 
(iodides) and some crystals based on fullerenes comprise two contributions at low temper- 
atures: a broadband ascribed to selftrapped excitons and, at higher energies, narrow-line 
(-band) optical transitions, fluorescence and stimulated emission with sharply structured 
absorption and emission transitions and a line-narrowing with inhomogeneous broadening 
due to free excitons. Spontaneous emission reflected in fluorescence or luminescence spectra 
was first observed as a free-exciton resonance luminescence coexisting with the broadband 
luminescence of the selftrapped excitons excited by electrons or synchrotron radiation in rare 



gas solids and alkali halides (iodides) [2lJ and later by other authors [22|, |23|. Transition 
from superfluorescence to amplified spontaneous emission with lost coherence was studied 
on super-oxide ion in potassium chloride in the intermediate regime as a function of temper- 
ature 23(]. Similar effects were observed also in structures based on Jahn- Teller (JT) active 
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molecules, e.g., fullerenes 6*60 24j . Laser-induced fluorescence line narrowing experiment 
of Cqq- hydro quinone crystals demonstrates inhomogeneous broadening of the zero-phonon 
lines. Similar observation has been done on Dicke superfluorescence in KC1 25| and the 
broadening was ascribed to unspecified internal strains. The peculiar high-energy narrow 
band luminescence in the exciton and Jahn- Teller bearing structures has been proposed to 
be caused by the emission from the localized "exotic" excited states within the phonon- 



selftrapped spectrum by Wagner et al. 
states proposed by Slonczewski [281] ). 



26 



271 ] (these states correspond to the tunneling 



Direct observation of the dynamic Jahn- Teller effect in Cgo 29( brought an evidence 
of important role of the ionic motion in the electron-phonon interaction. The quantum 
fluctuations led to formation of the tunneling states between two energy wells of the distorted 
phase (selftrapped electrons) which pertain to the undistorted geometry and recover free 



electron motion found theoretically [26|, [271, [3fJ . This picture corresponds to the concept of 
coexistence of selftrapped states and tunneling to localized (exotic) states at higher energies. 

Physical interest in JT models is motivated by importance of some spatially anisotropic 
complex structures (fullerenes, manganites, perovskites, etc) which exhibit JT structural 
phase transition. The transition is due to large local distortions related to the electron 
selftrapping. The JT model is a prototype model for electron coupling to two degenerate 
intramolecular vibrons of the same frequency (see ([6]) below) removing the degeneracy of 
electron levels (a-term) and phonon- mediated tunneling between the split levels (/?-term). 
Respective vibron spectrum is complex with level crossings and avoidings and needs exact 
diagonalization of matrices of confined level manifolds and proper statistical methods of 
evaluation. In crystals exhibiting high spatial anisotropy the rotation symmetry of JT 
molecules (a = (3, Q± = Q2, ^ being frequencies of the vibron modes) is generally broken. 
Therefore, it is reasonable to investigate the JT model assuming different coupling strengths 
for the onsite intralevel (a) and interlevel (/?) electron-phonon couplings and different vibron 
frequencies. Such a model can also be considered as a generalization of the exciton-phonon 
or dimer-phonon model, assuming the electron tunneling to be phonon-assisted. It implies a 
competition of the selftrapping distortion mechanism (localization) and a quantum tunneling 
(derealization) one. 

We will show that the lattice JT Hamiltonian of molecules with degenerate electron levels 
coupled with two non-equivalent boson (vibron) modes is formally equivalent to the Dicke 
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Hamiltonian generalized by adding a boson mode coupled with the level spacing. We shall 
investigate their excited spectra and related statistical properties vs those of the one-boson 
Dicke model. An impact of the additional boson mode which removes the double degeneracy 
of the JT electron levels close to the quantum phase transition from normal to supper-radiant 
phase reveals as a finite-size effect of level correlations mediated by the additional boson. 
We will study analytically an interplay of quantum fluctuations and symmetry breaking 
and numerically related characteristic signatures of a quantum chaos - the level spacing 
probability distributions, spectral entropies and spectral densities. We will demonstrate 
the localization-delocalization transition of the wave functions within the available lattice 
at the point of the phase transition between the quasi-normal and radiation phase. As for 
applications, we shall propose a zero temperature mechanism of the coherence dephasing and 
the broadening of the spectra, including the broadening of the zero-phonon lines. Further, 
we propose a similar impact on the quantum localization-delocalization phase transition 
from the Mott insulator to the superfluid state of ultracold atoms in optical lattices. 

Let us first shortly summarize the basic attributes of the Dicke model in order to motivate 
our use of the collective Holstein-Primakoff bosonization method for the two-boson JT model. 
The exciton-boson (photon, phonon) model of a lattice chain of N atoms is represented by 
the Hamiltonian 



where a\ n are spin variables; n is a chain site, N is the size of the chain, k = ±nl/N, I = 
0,1, . . . , N is the wave vector of the phonon wave, Qfc is the photon (phonon) frequency 
of the k-th mode, uo is the level spacing and A is the level-phonon interaction strength. 
Phonon creation and annihilation operators a^, satisfy conventional commutation rules 



allows also for the correct account of the transition to the thermodynamic limit iV — > oo. 
The Dicke model is obtained from the model ([1]) when the boson is supposed homogeneous 
(of long wave length) within a finite set of lattice sites, \k\ <C vr. Then, the collective spin 




k n=l v n=l k 



[ak,a~^] = 1. The coefficient l/\fN arises from the dipole coupling strength 7]. This factor 




over all atoms in the chain can be 



H = Via) a + uqJz + A^=(a f + a)(J+ + J J) . 




(2) 
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The collective (N-dimensional) spin variables J z , J + , J_ satisfy the same commutation 
relations as the individual spins: 

[Jz,J±] = ±J±; [J + ,J-] = 2J Z (3) 

1 i 

with J x = ~{J+ + J-), J y = —-^{J+ — J-)- It is handy to represent the realization of the 
spin operators of higher dimensions by the Dicke states. Each Dicke state is labeled as \ j, m) 
with the discrete index m ranging (at fixed j) from — j, + to j. The spin operators 
in the Dicke space are defined as follows: 



■Jz\j,m) = m\j,m) ; J±\j,m) = y/j(j + 1) - m(m ± l)|?,m±l) ; J 2 \j,m) = j(j+l)\j,m) 

(4) 

The subspaces of the Dicke states labeled by j are independent and, hence, one can consider 
them separately keeping j fixed. For the chain containing N atoms possible values of j range 
as 0, 1, . . . , N/2 for N even and 1/2, 3/2, . . . , N/2 for N odd; for given j the operators Jj are 
matrices with dimensions 2j + 1 (representations of the 577(2) group). The parity operator 

II = exp{i7r(a^a + J z + j)}, (5) 

commutes with Hamiltonian (j2J), i.e. the number of excitation quanta (a^a) + J z + j is 
conserved within a phase. Namely, the space of states is separated onto two subspaces with 
the number of excitation quanta which is either even or odd 0, Q] • Then the split subspaces 
can be considered separately. This conservation is broken at the quantum phase transition. 



33j, 



34 ] 



Two-level pseudospin-vibron (photon) JT model Hamiltonian investigated in 32 
has the form of the 2x2 matrix which in the pseudospin notation is written as 



H = Q(a\ai + a 2 ^a 2 + 1)1 + a(a\ + ax)a z — j3(a 2 ^ + a 2 )a x (6) 

Hamiltonian ([6]) describes the electron-boson interaction in a two-level JT molecule with 
two boson modes a, of equal frequencies Q (I in the first term is the unity matrix). In (jBj), 



the a-coupled mode lifts the double degeneracy of the 



diates the electron tunneling between the split levels 3J|. Different electron-boson coupling 



ocal JT level and /5-coupled mode me- 



strengths a and (3 present a generalization with reference to the conventional (rotational 



symmetric) JT model where a = f3 34]. 



In a similar way as for the exciton model ([T]) a natural generalization of the two-level 
pseudospin-vibron (photon) JT Hamiltonian can be done by extending the dimensionality of 



the (pseudo)spin space. Namely, the generalized model considered within the present work 
is characterized by the following Hamiltonian: 



H = fii Q + a\a^ I+n 2 Q + 4a 2 ^ J+-^(a t 1 +a 1 )-2J, + -^(4 + a 2 ) (J+ + J_ 



(7) 



The operators J are defined in representation of the SU{2) group Q in a 2j + 1- 
dimensional space spanned by Dicke vectors \j, m) as described above (Eq.(jlj)). We assume 
a general case of the rotation symmetry breaking with different boson frequencies Q± and 

In the molecular JT model the bosons represent the intramolecular vibrons. The rota- 
tional symmetry is broken by taking different frequencies and interaction parameters a ^ (3. 
The generalization acquires a twofold interpretation: (i) The JT lattice model, j lattice sites 
represent a finite number of identical JT molecules, i.e. two-level molecules coupled to two 
intramolecular vibrons (of the same intramolecular structure). Then the electron collective 
coordinate can be introduced representing a collective spin of iV two-level molecules (en- 
semble of qubits) coupled with two local vibron modes identical for each on-site molecule, 
(ii) The excited on-site molecule with two intra-molecular vibrons can be considered by tak- 
ing a set of many excited (j) levels of the molecule. (|7j) can be then interpreted as usual 
many-level system (a single atom with a spin > 1/2) coupled to two boson modes. 

The above concept of JT model can be applied, e. g., for the formulation of the JT 
polaron model in nanostructures where due to their finite size the quantum effects are 
rather detectable for observation. 

Additionally the photonic JT model 10, lUl] with collective boson modes homogeneous 
within the chain of atoms is within the scope of (j7|) as well. 

In the Dicke model, two bosons represent optical fields of the wave length larger than the 
size of the lattice and so they enable to introduce the collective variables for the electron 
levels. Mathematically, both the generalized Dicke and JT models are then equivalent and 
represented by equation (7). 

Thus an impact of the coupling of the level spacing with the additional boson a\ compared 
to the conventional one-boson Dicke model (j2J) is to be studied. 

Likewise in the exciton model there exists the same operator of parity II = exp{i7r(c4a2 + 
J z +j)}, which commutes with Hamiltonian (j7]). For the JT molecule with only two electron 
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levels (j 
p = ±l 



31 



J2) and equal boson frequencies this fact reflects itself in the conserved parity 



321 ] as an additional good quantum number. For parities +1 and —1 the spectra 



of that model are identical. As we shall see, for the present generalization this is the case of 
even number of molecules, that is for the main Dicke numbers j = 1/2, 3/2, .... 

Finally, a diagonal term of the form A J z can be added to the (j7j) where the parameter A 



accounts for the intrinsic level separation. This term accounts for a magnetic field [11]. In 
the JT model, this term is known to affect the shape of the adiabatic potential: it lifts the 
conical intersection of upper and lower adiabatic surfaces. 

Number and relative values of the pseudospin-boson couplings involved are of essential 
importance for specification of statistical properties of the pseudospin-boson models: while 
the excited spectra of various pseudospin-one-boson models were found to exhibit Wigner 
chaos when approaching the semiclassical limit , the excited spectra of the 

molecular JT model (Oj = Q o) ex hibit nonuniversal chaotic behaviour between the Poisson 



and semi-Poisson limit 
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35(| . The antisymmetric against reflection level splitting mode 



a\ of the JT model increases by one the number of local degrees of freedom when compared 
to the exciton model. Simultaneously, this phonon (photon) mode mediates additional 
level correlations and quantum fluctuations and thus it weakens the level repulsions. The 
interaction effects are modified by the frequency difference and are expected to imply an 
impact on the statistical properties of the lattice models such as quantum phase transition 
to the super-radiant phase and proper statistical characteristics of the quantum chaos. 

In Section II we investigate numerically and analytically the ground state and related 
phase transition of the extended Dicke model or, equivalently, of the lattice JT model (JTj) 
with two bosons of different frequencies. For analytical calculations we use the Holstein- 
Primakoff bosonization Ansatz for the pseudospin variables. We analyze the macroscopic 
phases for three involved oscillators related to the competing correlations (self-trapping and 
squeezing) due to the additional mode. In Section III we present analytical investigations of 
the interplay of quantum fluctuations and the symmetry breaking as a useful insight into the 
scenario leading to the macroscopic phase transition in the system of three oscillators. For 
a finite lattice, within the quasi-normal domain a medium quantum domain is shown to be 
created in which both the "normal" and "radiation" phase coexist. In Section IV we present 
numerical results on wave functions and statistical properties of excited states, namely level 
spacing probability distributions for various number of lattice sites, values of frequencies and 
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coupling parameters and compare them with those of the one-boson Dicke model. Further, 
the distribution of occupation of excited states over the lattice sites is illustrated by the 
state entropies as functions of quantum numbers of even states. The spectral densities as 
functions of energy are calculated as illustration of the effect of the interaction with the 
boson modes on the spectral density on the Dicke spectral space. The wave functions and 
the statistical characteristics reflect the medium quantum domain as being mixed from the 
localized and delocalized contributions of the normal and radiation phase, respectively. 

For presentation of some numerical results we use the resonance case tt% = Q 2 except for 
the cases where the difference of the frequencies makes qualitative changes. Throughout all 
the paper we shall use the term "radiation" phase for our model with suppressed coherence 
instead of the term " super-radiant" standardly used for the strictly coherent one-boson Dicke 
model. 



II. GROUND STATE AND PHASE TRANSITION 



Use of the method of collective pseudospin variables for two-level electron-boson JT and 
Dicke lattice models is appropriate under the conditions presented in Introduction. The 
collective coordinates enable solving these nonintegrable models efficiently, though approxi- 
mately. In what follows, we apply them for the investigation of a lattice of JT molecules or 
the two-boson Dicke model. 

Hamiltonian of the system can be written as a generalized version of ([1]) (k = 0) as follows 

H = Vl 1 Q + aJa^ I+Q 2 Q + 4«2^) / + ^-^j[«(a t 1 + ai)24 4) + /?(4 + a 2 )(4 ) + s ( ! ) )] , 

(8) 

where j assumes values in different subspaces j = 0, 1, 2, . . . N/2 for N even and j = 
1/2, 3/2, . . . N/2 for N odd. In what follows we will always take into consideration only 
the subspace with largest j = N/2 for each number of atoms in a chain. The frequencies 
of both modes are Q\, Q 2 and corresponding coupling constants are a, (3. Introducing the 
collective coordinates Jj defined by into equation (jSj) in the same fashion as for jTJ we 
arrive at Hamiltonian (j7]) representing either the generalized Dicke or JT lattice model. 
For j = 1/2, the local (molecular) JT model (with broken rotational symmetry if a ^ (3, 
Qi 7^ Q2) is recovered. For completeness, the standard Dicke-like model is the one-boson 



S 



version of Hamiltonian (|7j) without assistance of the bosons 1. 

For analytical treatment it is convenient to bosonize the spin variables in the two-boson- 
spin system ([7]) using the Holstein-Primakoff Ansatz 



j z = tfb - j, J + = tf^/2j-tfb, J- = y/2j-Vbb (9) 

for the collective pseudospin operators. Here, b are fictitious boson operators, [b,b>] = 1. 
This representation of the spin algebra preserves exactly the commutation relations ([3]) and 
makes it possible to convert the system to two or three coupled quantum oscillators as will 
be shown below. 

Linear terms in the Hamiltonian ([7]) can be excluded by the use of the coherent state 
representation with boson displacements chosen variationally. Let us displace the operators 
involved as follows 

a\ = c\ + y/ai, a\ = c\ + y/a2, = d) — \f~5~ . 



(10) 

Applying the Holstein-Primakoff Ansatz (Q, setting (flOj) into (J7J) and expanding the 
square root expressions in (Q to first order terms in Sd yields the form 

H = Qi ( i + c\c^j + Q 2 Q + c|c 2 J + ftiax + Q 2 a 2 + Qi^/ai (c\ + c^j + f2 2 \/o2 (c\ + c 2 
c\ + ci)d ] d + 2^fa[S d - V5(c\ + ci)(c/ f + d) - 2^/a l ~6(d^ + d) + (5 - j)(c\ + c x ) 



+ 



2a 



2j L 

+ 2^1(5-3) 
1 



(4 + c 2 + 2^a~ 2 ) • k \ -2^5 + (1 - 5/k 2 )(S + d) + ^-dU 



2k 2 



-(d ]2 d + d ] d 2 ) + V~5{{d) + d) 2 - 1) L k = ^/2j^~5 



(11) 



From the condition of elimination of the terms linear in boson operators from (fTTj) three 
identities for the parameters a« and 5 follow: 

2a 



fiiv/aii 



(s-j) 



j3y/5 4j -26 + 1/2 
72j VW^T~ 



(12) 



y/S(5-j) [4tt 2 a 2 -2tt 1 p- 



A] -25+ 1/2 



0. 
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It we consider a finite lattice (J finite), this set of equations implies three solutions 
provided Q 2 o? ^ flj /3 2 (1 + 1/8 j) as follows: 

1. v^= — v 7 ^, a 2 = 0, 6 = , (13) 

2. a 1 = 0, ^=£y/v(l + l-Y 6 = j, (14) 



^2 v V 4j 

3- V5T=-^-V^7(l-2^), v / ^=^^V / ^(l-^)> ^ = 2j(l-/i), (15) 

^2 ^ i ^2 

where in the last set j2 = — r- = — -, =-r- < 1, and a > f3 = \ZQ,i/Q, 2 P- 

8j(a 2 Q 2 -/3 2 1M 8j(a 2 -/? 2 ) 

As will be seen in the next subsections, first two solutions in the thermodynamic limit 
j — > 00 account for the common "normal" and "super-radiant" phases of a system, the 
transition between them at ai/H^ = describing the quantum phase transition (QPT) 

similar to what was reported for the Dicke models of superradiance jf], Q, 0, The third 
solution in the thermodynamic limit becomes trivially degenerated and valid only on the 
line of the said QPT. Nontrivial consequences for this solution arise for finite j only: From 
Eqs. ( TT3|) and (fl4"|) a transition to the rotational symmetry of the two-dimensional vibron 
oscillator is evident at a 2 = oi\ at the critical point a = /3\/tt-(1 + l/4j). In the limit 

V "2 

of large j — > 00 this is the point of the phase transition described below. This limit of 
higher symmetry exhibits qualitatively different properties when compared with the case 
a 7^ (3 \J Vli / Vl 2 (l + 1 / 4 j ) and j finite as we shall see in Section IV. Let us note, that the 

factors like ~ 1/j are the origin of specific quantum effects such as the squeezing of the 

4/3 2 

frequency lu 2 = — — cosh(4r), cosh4r = 1 + l/4j, of the oscillator (polaron) 2 (subsection 

li 2 

B). The effect of squeezing gets implications for the level spacing probability distributions 
especially in the case 3 when all three oscillators are coupled (Subsection C). 



A. Case 1: Polaron self-trapping in the normal phase. 



The solution ([13]) represents the normal phase with average zero number of macroscopi- 
cally excited electron level bosons 5 = 0. Hamiltonian (TTTT) related to this solution reads 

Hi = fix (~ + eld) + 2 Q + e|c 2 ) + ^(dtd - j/ 2 ) + /3(c 2 t + c 2 ){<# + d) 

+ ^(4 + c x )d)d - — (4 + c 2 )(d) 2 d + d)d 2 ) . (16) 
V 2 3 4j 
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Here, the level oscillator acquires the frequency of a polaron- dressed oscillator u\ = 
4a 2 /Qi which results from the self-trapping by the additional mode in contrast with re- 
spective factor of the Dicke model uq. Explicit interaction of the oscillator modes with the 
level bosons contribute only to excited states oc 1/y/j in ffTB]) . Thus, in the limit j —>■ oo 
Hamiltonian (|T6l) recovers the normal phase Hamiltonian of polaron level oscillators. For 
large j only linear terms persist (first line in (fTBT) ) and it can be easily diagonalized by a 



rotation in the plane of coupled operators c 2 and d 



71 



to yield the form 



H ld = n lC \ Cl + e 1;+ ClC 2 + e^ClCs - ^-J (17) 



with new effective oscillators C 2 , C3. The excitation energies of the system e\^± are 

Tl/2. 



^ ± 2 



(Q 2 2 -u; 2 ) 2 + 64a 2 ^ 



(18) 



From (|T8l) . the solution for e\ exists provided lo\ > 4j3 2 /Q 2 = u> 2 , or if ay^ > /?v^i- 
This phase is identified as the normal phase of the Dicke-like model without macroscopic 
excitations, Fig. [TJ The effect of different initial frequencies Qx 7^ Q 2 is just the shift of the 
phase transition point supporting the phase with the smaller of fij. 



From (fT6l) . critical Hamiltonian Hf u = Hi{(3 c = a-^/^M) at the point of the phase 
transition yields 



Hf* = Jni + (c*C a + 1/2) + ^(cjd + 1/2) - ^ . (19) 

The coupled undisplaced oscillators a 2 and b form an effective single oscillator of the 



frequency + {Aa 2 /^). The mixing results from the coupling due to the assisting 
oscillator a 2 between the levels split by the oscillator a\. The additional coherent oscillator 
ci remains decoupled. 

B. Case 2: Squeezing in the radiation phase. 



The second solution (J14j) identified further with the superradiant phase, can be treated 
on the same footing. For this solution the level bosons are displaced by 5 = j that is 
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they acquire a macroscopic number of excited quanta. Respective Hamiltonian (TIT]) yields 
1 + \ (\ 



H 2 = fii ( - + r\r. 



{ -h { ^ + c\c 2 ) + 2 ' 



2 ' \ - 

n 2 v + 4j 



Sd 



->/2a(4 + ci)(d t + d) 



+ ^ ( 1 + (rf f + df + 4^(4 + c x )d)d + -4=(4 + c^d + 
^2 V 4jy v/2j V2j 







(4 + C 2 ; 



V<7 



(20) 



Linear part of f)20p relevant in approaching the thermodynamic limit j — > oo represents 
the level polaron coupled with the oscillator 1. Quantum fluctuations squeeze its frequency 
by the squeezing parameter cosh 4r = (1 + 1 / 4 j ) . Transformation of this part of fl20l) by the 
unitary operator of squeezing S = exp[r(<fl" 2 — d 2 )} using the identities 



dtd = SdUS' 1 = d f dcosh4r + sinh 2 2r + (d 2 + d t2 )- sinh4r, 



•2r 



(21) 



yields the squeezed polaronic oscillator with renormalized frequency u 2 = 4/3 2 /Q 2 cosh(4r) 
and interaction k = \/2ae 2r , 

(\ \ ~ ~ 23 2 

-tt 2 ( - + c\c 2 J +cj 2 d t d-fi;(c t 1 +ci)(d t + J)--^-j-t ( ;2sinh 2 2r. (22) 



#2 = ^1 I g + r i r 

Diagonalization of (122]) implies three independent oscillators Ci, C3 and c 2 , last one remain- 
ing free, 



H 2d = e 2 , + C 3 t C 3 + e 2) _C 1 t C 1 + fi 2C Jc 2 - 



2/3 2 
CI 



-J 



r.:, = - 1 ttl + u 2 2 ± [(fi 2 - c^ 2 ) 2 + 64a 2 /5 2 e 4f ' fil1 ' 



,2 

-2,± 



where 

^ 1 izj + ^ ± |^zj - c^f + b4ar/re" ^ 

Energy of the squeezed ground state is H 2 q = —2(3 2 /VL 2 ■ j — a; 2 sinh 2 2r. From 
obtains the condition for stability of the radiation phase 



of 



£->-< 

^2 ^1 



(23) 

(24) 
one 

(25) 



Let us remind that for j — > 00 the squeezing parameter r — > (sinh4r ~ 1/^/2]). 

Close to the phase transition point (3 C = a^Q 2 /Q\ where e[ ^ = 0, we can find the energy 
from the side of the normal phase a — > /3 C _ 



=(-) 



. a: 



3/ 2 fif 



V2 (n 2 n 2 2 + 16a 4 



1/2 



1/2 



(26) 
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From the side of the radiation phase, at (3 C+ <— a+, the energy reads 

c(+) , s^gV- (Q f7T 



y/2 (ttjttl + 16pA 



1/2 



(^-av^) • (27) 



In the limit j ; — >• oo the squeezing parameter r vanishes, and the energy dependence becomes 
symmetric with simultaneous interchange a <-> /3, fii <-> f^. Note that in the resonance case 
Oi = Q 2 the cusps in ( 1261) . (I27j) are symmetric around the critical point a = (3 (Fig. [1]). 
However, for finite j this symmetry is broken because of the factor exp(2r) in the numerator 
of ([27)1 : the radiation phase (1271) is suppressed by the squeezing when compared to the 
normal phase. The non-symmetry of the branches below and above the critical point is 
also evident from numerical results for the excitation energy from the exact Hamiltonian (J7J) 
presented in Fig. [TJ Let us note, that the critical exponent 1/2 near the critical point in (J26J) 
and (1271) is characteristic for the mean field second order phase transition. The branches 
(|26|) and (f2"7|) represent a generalized version of the branches by Emary and Brandes Q : The 
analogy becomes obvious if we set A c = \^/ujuJq = J t^Pc, with ui = fij = fi 2 , cu = 4ct 2 /f2 1 . 

V \l 2 

In the limit j —>■ oo, similar critical behavior as (Tl9|) can be found for the energy of the 
phase 2, H^ rit = H 2 U h^-f3 c = a). From (EHJ, one obtains 

V "2 

H? u = \^M^\c\C x + 1/2) + 2 (4c 2 + 1/2) - ^j, (28) 

where now the displaced oscillator c 2 remains decoupled and the undisplaced oscillator C\ 
mediates a mixing with the split level polaronic oscillators of the frequency uj 2 = 4j3 2 /Q 2 . 

The analysis of both phases shows that within linear approximation the oscillator c\ in the 
radiation phase plays qualitatively the same role as it does the oscillator c 2 with simultaneous 
interchange of the polaron frequencies LO\ <-> uj 2 (coupling constants a <-> fi\J% 2 - Hence, in 
linear approximation two oscillators mix to two-dimensional effective oscillator while the 
remaining one is decoupled. For the resonant case Q± = Q 2 and a = (3 {uj\ = u 2 ) a 
spontaneous transition occurs (in the plane of vibron coordinates (Qi,Q 2 )) to the rotation 
symmetric two-dimensional oscillator of a qualitatively different behavior which is out of the 
scope of this paper. This case implies an abrupt change of the statistical characteristics of 
the quantum excited levels (section IV) related to the higher symmetry of the problem 33]. 
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C. Case 3. Intermediate domain of mixed quasi-normal and "radiation" phase 



The solution (fi~5l) . valid for finite j, implies for Hamiltonian (fTTj) following form 

2a^l - p,(c{ + ci)(d) + d) + (5 ^~ (c| + c 2 )(d f + d)+ 

^(i - p){d) + d) 2 + -^(cj + Cl )dU + ^ ^2 (c t + cs,)* d 
+ J^ VT z Ji (ct + c 2 )(j + dy 

2V2J v 7 / 1 

_ -^=(i - fiy/^d + dt^) - -A=(cJ + C2 )(d t2 d + d f d 2 )- 

v (1 - 2/i) + -^r M1 - a — ^ ^ — ' (29) 

where // = — p r- < 1, a > p\ — . 

8j(a 2 n 2 - ^nA V sfe 

From (T291 one obtains for the quasi-classical oscillators (c| + ci) = qiy/2£li, (c 2 + c 2 ) = 
q 2 \p2£hi and for the level polaron (d* + d) = Q^8a 2 ft/ili dynamic equations 



a 



q\ = -n 2 qi + 4a Jl - p, y/u^Q - -=w V n iQ . (30) 

g 2 = -i2 2 g 2 - 2/3 — V^i^Q r - t= <2 , (31) 

Q = -cu 2 <5 + 4av/l - ^A/w^igi - 2(3^—— -^/uAl 2 q 2 - -^— ^ J 1 uQ 

~ P\ — —u\f£l 2 q 2 Q — % ^ _f u)\[Sl 2 q 2 Q 
V 3V v\7 ft 

1 ^ 'l-fi^Q+W^rt,. (32) 



Linear approach used for the normal and radiation phase far from the critical region can 
not be applied for the mixed critical region; Instead, nonlinear dynamic equations are to 
be solved. The dynamics of the level oscillator Q in (I3"2"|) is influenced by the squeezing 
parameter ft: If j is sufficiently large, ft is small (a 2 > j3 2 fli/Q 2 ) and uj = 4ct 2 ft/fii <C Q 1 . 
It follows, that we can neglect adiabatically the intrinsic dynamics of the oscillators q\ and 
q 2 and suppose it implicitly determined by the dynamics of the squeezed level polaron. 
(We confirmed numerically the instability of the oscillations of qi and q 2 ). Hence, we can 
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eliminate the oscillators q% and q 2 from (|32|) and for the squeezed level polaron Q, in the 
limit j — > oo, i.e. — > 0, we obtain the dynamic nonlinear equation for the classical "order 
parameter" , 



16a 4 / P 2 \^ 3 • 2 4 a 4 ^ v / 2(l - P 2 /a 2 )^ 2 2 6 a 6 (l - J?ja 2 ) , 



° = + V 2 Q ~n^ Q ' (33) 

or, using the transformation Q = q H 1 the quadratic term is eliminated 

2a V2(l-/3 2 /« 2 ) 

and one obtains the normal form 

q = Aq- Bq 3 + F, (34) 

where 



„ _ 8a 2 (2a 2 + / 3 2 ) D _2 6 « 4 (« 2 -^ 2 ) n 40^/2(1-^) 8 ^3 
A = —k , B = — , > U, -r = -— h 



(35) 

Contrary to the previous cases, the transition described by equation ( 1341) between the normal 

and radiation phase is of the first order and supports a coexistence of both phases in the 

sector of the normal phase a > (3, or u\ > uj 2 . We come to the conclusion that in contrast 

with the Dicke model the presence of the additional oscillator opens a sector of a mixed 

phase with partial occupation of the excitation space of all three coupled oscillators. A 

B A 

transition between two non-equivalent minima of the potential V = — g 4 — — q 2 — Fq arises 
and can be considered as the oscillation-assisted tunneling (hopping). The formation of the 
radiation phase is driven by the force F which is compensated by the nonlinearity so that 
the "nuclei" (bubbles) of the radiation phase stabilize for the set of parameters (|35|) . The 
exact solution to ( l34l) is given by 

n x + cos(wr) 

Q = a — 7 y (36) 

n 2 + cos{wT) 

where 



A nl _ . nl - 1 F 



n 1 = (2-n 2 )/n 2 , a 2 = ^jfrf' ^ = 2A ^+2' a = ~^ 2 + n2) - (37) 



The solution ( 1361) represents nonlinear non-sinusoidal periodic oscillations. From equation 
(|33|) it is evident, that the nucleation process vanishes identically when approaching the 
transition point ui\ — ui 2 , where the effective polaron localization energies W\ and u) 2 in both 
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sectors are equal. Let us remind that the scenario described above can be perturbed by the 
fluctuations neglected by the adiabatic elimination of other two oscillators (qi and q 2 ) and 
by excitations due to the finite-size effects. 

Let us note that, in spite of independence of Eq. (1551) on j in the thermodynamic limit 
j — > oo, the dynamic variable Q is j— dependent so that the resulting mixed phase exists 
only for finite j as expected. 

Numerical evaluation of the order parameter J z from the exact Hamiltonian (J7J) presented 
in Fig. [2] illustrates the existence of the intermediate phase by non-zero occupation in the 
neighborhood of the point a = (3. 



III. INTERPLAY OF QUANTUM FLUCTUATIONS AND SYMMETRY BREAK- 
ING 

A qualitatively new situation occurs in the nonlinear regime at finite j, when all three 
oscillators couple via interplay of fluctuations due to finite j and of breaking the rotational 
symmetry (\a — (3\ > 0). 

As the next step, we shall include in (fill) small terms up to 0(1/ j) to find the effect of 
quantum fluctuations. To elucidate the role of fluctuations and the difference between one- 
and two-boson cases let us apply first the method we will use in what follows to the Dicke 
model (j2J). We displace again the operators (see 

o+ = c f + Va, 6 f = S - V5 (38) 

to exclude linear terms from Hamiltonian (j2j). Here, a and 5 are the mean values of excitation 
number (a T a) and (b'b) of related bosons a and b, respectively. Putting the Ansatz (I3"8"j) into 
(J2J) and eliminating linear terms, one obtains two known solutions [3] 

1. a = 0, 5 = 0, normal phase, "1"; (39) 

2\ l l 2 / n\ V 2 



2- ^ = ^ ( 1 - {0) j , V5 = (l - ^) , radiation phase, "2" (40) 

which determine displacements of the normal phase "1" with zero number of macroscopically 
excited bosons, and of the radiation phase "2" at A 2 > with both numbers ~ j. Then, 

Hamiltonian (j5J) for the normal phase is reproduced, 

Hid = ^(c f c + 1/2) + u Q {rfd - j) + A(c f + c) U + d J J . (41) 
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Let us write Heisenberg equations for the phase 1 

S 2 d + Sd 2 



it = nc + A \d) + d ^ j , (42) 

id = co d + A(c f + c) - ^(c f + c)(2dU + d 2 ) . (43) 

Within the adiabatic approximation, we can neglect the intrinsic dynamics of the level 
oscillator d, assuming plausibly Q <C ujq. Thus, we suppose d(d) — d)/dt = in f j43l) so 
that the oscillator can be eliminated from this equation and put into (142 p . Defining then 

operators q = (c^ + c) and ft = i\ I — {c* — c), the dynamic equation for the coordinate 

v 20 V 2 

q = (q) up to the order 1/j can be written: 

q = -n 2 ( 1 - %) q - + 0(1/ j 2 ) . (44) 

4A 2 . ... . . . „ „ / 4A 2 ^ 1/2 



If — — < 1, there result oscillations with the squeezed frequency dp = CI [ 1 — 
about the displaced center. The quasiclassical potential corresponding to equation (jiljl 



Q 2 A 4A 2 \ 2 A 4 fi 2 lS 



yields 

V = V ( 1 " ^ ) ^ + 7^3 9 4 + Vo- (45) 
For each j, the potential (145p implies a second order phase transition at A 2 = ujqQ/4 
which has been found as the critical point of the phase transition between the normal and 
super-radiant phase in the Dicke model [t]]. For A < A c , the normal phase is recovered being 
stable while at A > A c the radiation phase is stable. Let us note, that the system (145]) 
is globally stable for each finite j and the phase transition is supported by the quantum 
fluctuations ~ 1/j. The method we present here is so verified as being able to reproduce 
basic results known from the standard approach used before |tJ and applied for the present 
model in previous Section. 

In what follows we shall perform analogous calculations for the system with two boson 
modes (JTBj) . For this case, Heisenberg equations related to ([IB]) for the phase 1 read as 
follows 

ic x = n lCl + ^LdU, ic 2 = n 2 c 2 + p(d ] + d)- —(d) 2 d + dU 2 ), 

V2j 4j 
id = ^-d + ^={c\ + Cl )d + (3(c\ + c 2 ) - A (4 + c 2 )(2dU + d 2 ) . (46) 
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4a 2 

For strong interaction, it is plausible to assume io\ = — — ^> f^, Q 2 and, consequently, 
to apply the adiabatic approximation neglecting intrinsic dynamics of the level polaron d 
of high frequency lo\, d — 0. Consequently, it can be eliminated but it determines the 
dynamics of the slaved modes qi = (2VLi)~ 1 / 2 (c\ + Ci), q 2 = (2Q 2 )~ 1 ^ 2 ( c 2 + °2) implicitely. 
Hence, we use the stationary expression for q = (8a 2 /n 1 )- 1 ' 2 {d^ + d) from flU and put it 
into remaining equations for q± and q 2 . Up to the lowest order term in one obtains 

* = — Aawr q2 (47) 

q 2 = -n 2 2 (l- ^— i g 2 - - gi? 2 - (48) 



We received the squeezed frequency of the oscillator q 2 which justifies again the use of the 
slaving principle for the oscillator qi, Qf ^> fi 2 (l — (3 2 /a 2 ). Its dynamics is then implicitly 
ordered by the dynamics of the oscillator q^. Then, similarly as in the previous case, we 
receive 

dq 2 

where 

fl 2 ( p 2 nA 2 P^jQ 



Equation fl50|) implies a continuous phase transition at a — (3 = ^^TjQ^ to the 
order 1/j for each j between the normal oscillators ( T4H1) with squeezed frequency CJ2 = 
Q 2 (l — /3 2 /a 2 ) , (a > /5), and a nonlinear sector analogous to the super-radiant phase of 
the Dicke model (ft > a). Note that the potential for the Dicke model (T4"5]) and fl50l) become 
identical if we set ujq = 4a 2 /f2 in equation fT4"5l) . 

We conclude that there occurs a sequence of phase transitions for each finite j between 
the quasi-normal and radiation phase (J3 > a) because of softening (squeezing) the frequency 
of the oscillator uj 2 = by the parameter (3 /a. 

According to Eq. ( 1261) and (1271) the critical energies at both sides of the phase transitions 
are symmetric when changing a «-> f3 and Qi ^ Q 2 . In the radiation phase, for j finite, the 
nonlinearity bears one-instanton solution 

</2i(r - 7-,,) : ±^i-./r£ ( 1 - ) tanh I % fl - £ V ' (r - r ) I (51 ) 



V P 
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Here, cfe = 92 y§^, ft = ft\J~^,- The one-instanton solution (15TT) is associated with the 
tunneling between the extrema of the potential inverted to fl50|) if ft > a and r — > it. Hence, 
at finite j's there appears a new instanton phase as a precursor of the phase transition at 
the maximum softening of the frequency at ft — > a in the radiation phase. More generally, 
there occurs a sequence of repeated tunnelings (oscillations between two equivalent minima 
of a local potential) for each lattice site. The same effect occurs in the Dicke model if 
4X 2 /(ujqQ) > 1 (T4"5l) . Moreover, in the present model, there exists the coupling between 
the oscillators 1 and 2 (I47II48P which was neglected in the adiabatic approximation (1491) . In 
fact, within more subtle calculations there would occur tunnelings mediated by two coupled 
oscillators (one of them being a polaron) for each j instead of one of the adiabatic treatment. 

The time of the tunneling T = cu^ 1 is determined by the squeezing of the frequency 



ojt = ( 1 — — 2 J , i.e. the tunneling ceases at the transition point. The frequency 
u>t determines the curvature close to the minima of the bistable potential floD]) . This initial 
frequency enters the probability of the tunneling from the ground state of the oscillator 
close to the minimum during the time r when starting at r (e.g., 38]). According to 



Dekker 



38| the trajectory of the tunneling particle in time exhibits chaotic features if the 



initial frequency is squeezed. This implies that the quantum tunneling should be associated 
with the chaotic features inherent to finite lattice or excited spectra of the standard Dicke 
model (1451) or our model in the present approximation. However, there is self-evident a 
question about the possible mechanism of chaos production due to the squeezing in this 
models. Going from the present considerations, instead of the strict tunneling (l5Tj) . we 
propose a more realistic mechanism of hopping, i.e., the above mentioned tunneling between 
the squeezed minima assisted by the neglected (when going from ( |46l) to ( |49l) ) nonadiabatic 
fluctuations of the oscillators c\ and d coupled by the nonadiabatic nonlinear terms ~ 1/ ^/j 
and similarly of the oscillator d ( H3i) for the Dicke model. We remark in passing that the 
mechanism of hopping is known, e.g., from the transport theory of disordered systems at 
finite temperatures [391 ] . 
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IV. STATISTICAL CHARACTERISTICS OF EXCITED LEVELS AND WAVE 
FUNCTIONS 



The distribution of the nearest neighbor level spacings (NNS) of excited quantum levels is 



the standard testing point for investigation of the issues in quantum chaology J36lj . Recently 
we have reported results on NNS distribution of the molecular JT model ([6]) 33j which is 
the case j = 1/2, Q\ = Q 2 , « 7^ /? in the notation of the present article. It was shown 
that these distributions essentially deviate from the Wigner distribution typical for well 
developed quantum chaos. Our results show NNS distributions nonuniversal but far from 
the cases with higher symmetry (rotational a = (3 and linear a > a < /3) they tend to 
a form close to the semi-Poisson distribution P(S) = 4S"exp(— 2S) already encountered in 
some problems, in particular in the Anderson model close to the metal-insulator transition 

n 

[40J • The peculiarity of the Jahn- Teller system with equal frequencies was that the Wigner- 
Dyson level spacing distribution has never been reached, and the said semi-Poisson form of 
the statistics appeared to be the most "chaotic" one. 

In this section we present statistical characteristics (NNS distributions, entropy of level 
occupation, and spectral density of states) of the spectra for fix 7^ Q 2 , j > 1/2. We solved 
numerically the eigenvalue problem for the quantum Hamiltonian ([7]). The Hamiltonian 
was diagonalized with the basis of the boson Fock states for bosons 1 and 2. When taken 
Nx and N 2 Fock states for each of 2j + 1 electron levels, there is produced an inevitable 
truncation error, so that the results were checked against the convergence (with changing 
the numbers N\, N 2 ). Only about ~ 1100 lower states were used for calculating the statistics 
out of typically ~ (8 4- 10) x 10 3 . The obtained raw energy spectrum had to be treated by 
an unfolding procedure in order to ensure the homogeneity of the spectrum (constant local 
density of levels). Thereafter the statistical data were gathered in a standard fashion from 
small intervals in the space of parameters (a, (5). The calculations were performed essentially 
along the same lines as our previous calculations for the Jahn- Teller problem with j = 1/2 



33l . |35[ where additional details were given as to the convergence check, unfolding procedure 
and gathering statistics. 

The results for the level spacing distributions for different phonon frequencies Qi 7^ fl 2 
and j > 1/2 show rather vast variety ranging from Wigner to Poisson distributions as 
limiting cases but recovering also the semi-Poisson distribution in a rather wide range of 
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model parameters. With increasing j, general tendency consists in increasing the diapason of 
parameters where the NNS shows deviations from the Poisson distribution towards chaoticity 
up to the maximum degree of chaoticity given by the Wigner distribution. So, for the case 
j = 1/2 and equal frequencies the maximal degree of chaoticity being expressed by the 
semi-Poisson distribution in accord with our previous results. For moderate j (j = 7/2 and 
Qi ^ Q2 exemplified in|3]there appears a well-marked domain of maximal chaoticity (Wigner 
distribution). For example, for f^/^i = 2 the domain of Wigner distribution stretches for 
the values of parameters 1 < a, (3 < 3. 

We can conclude that the presented results essentially deviate from the Wigner statistics 
of levels imposed by the one-boson Dicke model Q]. For comparison, we performed the same 
calculations of NNS distributions for the standard one-boson Dicke model (Fig. H]). The 
distributions follow closely the Wigner form of the P(S) curve, but a good agreement is 
achieved only for high values of j. For the resonance case Qi = Q2 the quantum chaotic 



statistics (Wigner surmise 36]) changes towards a non-universal with respect to a, (3 but 
tending to the semi-Poisson intermediate statistics achieved asymptotically far from the cases 
of the special symmetry (rotation, linear). In contrast to one-boson models, this universality 
is observed even for small number of electron levels and seems not to change with raising j. 
Therefore the impact of the second boson mode is that it essentially changes the quantum 
statistics. For the non-resonance case Qi 7^ Q2 we observe a considerable suppression of 
chaos manifesting in the reduction of the "pure" domain of Wigner chaos in the space of 
the parameters a, (3. 

To visualize the wave functions of the excited levels of the system we considered the 
spreading of the wave functions over the electron levels and integrated out two boson degrees 
of freedom. Thus we found numerically the probability P± of the wave function (ftf 1 ^ to 
occupy a given electronic level i of the system. The full wavefunction depending on both 
vibron and electron level variables is calculated by the same scheme of the diagonalization 
of the quantum Hamiltonian in the representation of Fock states as described above. If 
(3 = 0, Hamiltonian equation is split onto N(= 2j + 1) independent equations labeled by 
i = 1, . . . , (2j + 1), so that the wavefunction for each given state n is localized on some 
electron level i. If (3 7^ 0, Hamiltonian matrix is no more diagonal (since J x is nondiagonal) , 
so that its components get interacting. The resulting eigenfunction is in general (2j + 1)- 
fold function occupying each electronic level (see Chapter 1 as for the discussion about 
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normal and super-radiant phase in the thermodynamic limit). Let Xi {Q11Q2) be the i-th 
component of the (2j + l)-dimensional vector of the eigensolution of the Hamiltonian matrix 
equation for n-th energy level in the "coordinate" representation {Qi oc b\ + bi, 1=1,2). Then 
the occupation probabilities of the i-th electronic level are p/ = J J \Xi^\ 2 dQidQz- In this 
representation we exemplify the wave functions in Fig. [5] for the levels n — 1,4 and the 
parameters around the point of QPT a = 2, (3 — 1.95, 2.05 and Q± = ^2- The abrupt 
change in the shape of wavefunctions when going from normal to super-radiant phase is 
easily perceivable not only for the ground state, but for lowest excited states as well (in this 
example for n=4). However, for higher excited states the wavefunctions generally spread 
over the whole available electron space, irrespective to the relation between a and (3. The 
last relation is also reflected in the statistical properties of levels: they are invariant with 
respect to the exchange a «-> (3, fl\ <-> f^; likewise it was already shown for the generalized 



JT model 



331 ] . In order to characterize quantitatively the behavior of the wavefunctions of 



excited states we introduce the entropy of level occupation 

N 



Sn = -'£ P i n)l0 Z P i n) > ( 52 ) 



t=l 

N 



where Yl Pi = 1- I n Fig. [6] for Qi = Q2 we plot the said entropies for first 1000 quantum 

i=i 

levels of a system with four electron levels (j = 3/2). One can see from these figures that 
in the most excited states the electronic levels are equally populated (with probabilities 
1/(2 j + 1), e.g., last Figure in Fig. (jSJ)), thus the said entropies tend to their limiting value 
log(2j + 1). For four levels the entropy yields the value S4 = 2 log 2 = 1.38 which is confirmed 
by Fig. [6]for a 7^ (3. However, among these extended states there emerge "localized" levels of 
lower values of the entropy which are characterized by a relative localization of wavefunction 
on smaller number of electron levels. Such levels form marked branches seen in Fi g. El which 
remind of the branches of "exotic" states in the generalized Jahn- Teller model 34j]. Most 
of the states show an intermediate degree of localization between both of the limits. These 
intermediate electronic states correspond to the mixed intermediate domain of the normal 
and radiation phase in the respective boson space. 

Another way to visualize the complex structure of the excited states and corresponding 
wave functions is the representation of the spectral density given by the imaginary part of 

42]]. It shows the characteristic frequencies of the final state of the 



the projected resolvent 
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evolution of a system starting from the projected Hilbert space and returning to it. The 
spectral density of states F(E) can be defined with respect to some initial state |^o) (the 
ground state, see below) and is connected with the return probability to this state in the 
course of the system evolution through the exact states $ n , eigenfunctions of the energies 
E n of the system at small e ~ 0: 

F(E) = lm(9 \(E -H- fe)" 1 !^) = e £ { ^fj[ £2 , (53) 

where |$ n ) and E n are the eigenfunctions and corresponding eigenvalues of the Hamiltonian 
H. The small parameter e fixes the rules of handling the poles of the Green function in 
the complex space. In Fig. (J7|) we show the examples of the spectral density for the values 
of parameters a, j3 below and above the critical line. As the initial state \^q) we took the 
ground state with a — /3 — 0, that is the product of the electronic Dicke state \j, —j) e i and 
the phonon state |0, 0) p h with zero number of both bosons 1 and 2. Thus the spectral density 
in Fig. (JTj) characterizes the evolution of the system with the interactions a, (3 switched on 
in the initial time. The peaks at each E n measure the overlap between the initial state \l/o 
and the eigenstate $ n for e ~ 0. 



V. CONCLUSION 



The numerical analysis of the generalized (two boson modes with different frequencies) 
Dicke model or the equivalent generalized lattice JT model with fully broken local rotation 
symmetry (different coupling constants a / (3 and different frequencies fl± ^ Q 2 ) brings 
following new results: 

(i) The critical point of the second order phase transition , (JT]), is modified to a = 
/?y§^ when compared with the critical point of the one-boson Dicke model a = (3. The 
renormalization of the couplings by the frequencies is due a reformulation of the critical 
point in terms of frequencies of the dressed oscillators as u\ — 002, 0J1 — |p , oj 2 = 

(ii) The electron level order parameter (J z )/j © shows the finite size effect for finite j 
close to the phase transition point which smoothes the phase transition to a rather broad 
at small j "mixed" region. The effect is slightly nonsymmetric around the critical point 
supporting the order parameter in the normal phase. 
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(iii) The level spacing probability distribution Qof the model under consideration bears 
a drastic change from the Wigner-Dyson P(S) ~ S exp(—S 2 ) distribution of the one-boson 
Dicke model ([111 by a non-universal transition up to close to the Poisson distribution P(S) ~ 
exp(— S). For equal frequencies, Q\ = Q2 an d different coupling constants, however, the 
distribution ranges non-universally from a close to the semi-Poisson P(S) = 4S'exp(— 2S) 
to the Poisson distribution. Evidently, the change to the range between the Wigner-Dyson 
and the semi-Poisson distribution is due to the difference in the frequencies. 

(iii) In the electron space, the phase transition reveals in a drastic change of the wave 
functions from those spread over small number of levels (sites) in the "normal" phase to the 
macroscopic distribution of the excitation over all the levels (lattice) in the radiation phase, 
(|SJ). This behavior shows up a non- universality with regard to the interaction parameters 
and inhomogeneity throughout the spectral levels. 

(iv) Entropies of occupation of the electronic levels (jSJ) illustrate the measure of localiza- 
tion of the states throughout the spectra (or the lattice). There is evident the asymptotic 
line of the delocalized states at S = 1.38 and the localized states with various extent of 
localization. Close to the line from below there is evident a relatively high density of weakly 
correlated states. The lowest levels show up branches of the localized states. Let us remind 
that the exceptional case of a = (3 belongs to a higher class of symmetry out of the concept 
of the present paper. 

(v) A measure of the localization of the spectral states as a function of energy is the 
spectral density of states (I7j). It measures the evolution of the states from the pure Dicke 
states (at = 0) in t — to the exact numerical states of the system at interactions switched 
on at t > 0. During the evolution of the system to the final state the system overlaps with 
the initial state at the poles of respective propagator F(E) fl53|) . The peaks of the function 
F(E) measure the overlap between the initial Dicke states \&o an d the exact eigenstates 
$ n (£) for e ~ 0. 

Beside the exact numerical analysis, use of the Holstein-Primakoff Ansatz made us possi- 
ble to treat analytically, though approximately, the present model. The analytical treatment 
brings a useful insight into the mechanisms behind of the properties illustrated by the numer- 
ical results. We have shown that an additional boson mode coupled with the electron level 
spacing dresses the normal-phase oscillator to a self-trapped one similarly as it was dressed 
the oscillator in the radiation phase. This mode-dressing modified by squeezing cooperates 
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with nonlinear terms at finite j: It modifies both the normal and the radiation phase close 
to the transition point a = P\ tt-- The intermediate region of the "nucleation" occurs as a 

V "2 

consequence of the competition of the self-localization due to the polaron dressing modified 
by squeezing of the level boson mode and of the boson-assisted tunneling. Both quasi-normal 
and radiation phases coexist within the sector of the quasi-normal phase a > (3. We have 
shown that the formation of the radiation phase is driven by the boson-assisted tunneling 
(hopping) and is compensated by the nonlinearity due to the selftrapping so that the "nu- 
clei" of the radiation phase get stabilized. So, in contrast with the one-boson Dicke model, 
for finite j the additional oscillator opens a sector of a "quantum mixed phase" with partial 
occupation of the excitation space of all three coupled oscillators. 

We have shown that in the radiation phase there occurs a sequence of repeated tunnel- 
ings (oscillations between two equivalent minima of a local potential) for each lattice site. 
The radiation phase can be then represented as an almost ideal instanton-anti-instanton 
gas phase (instanton lattice). However, there exists a weak coupling with two oscillations 
assisting the tunneling because of non-adiabatic fluctuations due to the finiteness of the 
lattice. They are necessarily involved by the model and moderate the polaron dressing of 
the level mode fl29l) . As a result, they perturb the ideal periodicity of the instanton lattice. 
In a variable extent the quantum fluctuations due to the level correlations persist within 
the whole radiation phase and represent the dephasing mechanism of thorough coherence of 
the phase. We have shown that the sequence of local " second-order phase transitions" and 
so the instanton-anti-instanton gas phase exists in the standard one-boson Dicke model as 
well, in contrast to the mixed phase which does not exist in this model. 

Numerical calculations show reflection of the described three phases in the statistical prop- 
erties. We have demonstrated a drastic qualitative change of the level spacing probability 
distributions from the universal Wigner of the one-boson Dicke model to the non-universal 
transition from the Wigner distribution of the localized states of the quasi-normal phase 
through the critical semi-Poisson distribution of the mixed domain to a close-to-Poisson 
distribution for the almost ordered radiation phase. This change is attributed to the weak- 
ening of the Wigner chaos caused by the tunneling mechanism mediated by the additional 
mode. The fluctuation effects in the "radiation" phase are respectively weakly j-dependent 
for small j. 

For finite lattice systems, where the localization-delocalization phase transitions occur, 
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we suggest possibility to apply the present results to interpret finite-size fluctuation effects. 
We propose the applications of the presented results for finite systems where fluctuations 
destroy quantum coherence and a second order phase transition is modified by the tun- 
neling mediated by the fluctuations. The destroying the coherence of the radiation phase 
by the boson selftrapping is a source of the dephasing mechanism in the radiation phase. 
The spectral phenomena as is the inhomogeneous broadening of the zero-phonon spectral 
lines and the broadening of spectra are due to the selftrapping in the quasi-normal phase. 
This interpretation can be applied, e.g., as a mechanism responsible for the inhomogeneous 
broadening of the zero-phonon lines of the Dicke super fluorescence in KC1 {25 1. 

The finite-size fluctuations enable the formation of a critical phase of mixed localized and 
delocalized phases. An example of such a system with analogous picture of quantum phase 
transition provides the ultracold Bose gas in an optical lattice where the phase coherence of 
Bose atoms is destroye d by the suppression of the tunneling due to the localization by the 



strong lattice potential 42j , 



43| | . From the competition of the tunneling and the localization 



due to the lattice there results the phase transition between the highly coherent superfluid 
ohase of Bose atoms and the localized Mott insulating phase of the atoms on the lattice 



4=11 ] . If some fluctuations, e.g., due to lattice vibrations, mediated the tunneling within the 



insulating phase, as well as perturbed the coherence of the coherent superfluid phase, the 
system of Bose atoms in the optical lattice would exhibit analogous properties as those of 
the two-boson Dicke model. 

We propose also the additional mode as a mean to reduce chaos in related spectral 
properties of optical devices based on the Dicke model. 
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e 



a=0 . 5 




FIG. 1: Phase transition between the normal and radiation phase in the generalized Dicke and 
J-T lattice model for a = 0.5, 2 and resonance case fii = = 1. The symbols show the numerical 
results for the excitation energy of the first excited state for different j, the solid line in each 
figure is the analytical result for QPT (|18p . (|24p valid for j — > 00 . The cusp-like behavior at the 
critical point uj\ = 4a 2 /VL\ = 4/3 2 /Q 2 , = 0J2 appears already for relatively small number of sites, 
e.g., for j = 7/2. For small j, the fluctuations smooth the cusp especially at weak couplings. The 
non-symmetry about the critical point by the reduction of the energy in the radiation phase for 
a = 0.5 is due to the effect of quantum fluctuations (squeezing) and vanishes with increasing j and 
a. For the case of different 7^ O2 the picture is qualitatively the same, there occurs a shift of 
the transition point inwards the phase with higher Qj. 
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FIG. 2: Order parameter (J z )/j for j = 1/2,3/2,5/2 and Sli = f^- The finite-size effect (j finite) 
of the mixed phase about the transition vanishes for j — > oo. 
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FIG. 3: Level spacing distributions for j = 7/2 and different (a, (5) (parameters are scaled to 
fii = 1). Upper row ((a,b)): 0,2 = 2 0i; bottom row ((c,d)): = 0.5 Q%. The curves show 

Wigner (solid), Poisson (dashed) and semiPoisson (stars) distributions. NNS distributions in b), 

r I 

c) are close to the semi- Poisson distribution P(S) = 4Sexp(— 2S) [33|]; histograms in a) and d) are 
almost perfect Wigner (P(S) ~ S'exp(— S 2 ) and Poisson P(S) = exp(— S) distributions. 
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FIG. 4: Level spacing distributions for one-boson Dicke model ([2|) (boson frequency f2 scaled to 
1). At large j = 7/2, 11/2 almost pure Wigner chaos P(S) ~ S'exp(— S 2 ) is revealed. The smaller 
j, the stronger are level fluctuations growing with decreasing A. 
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FIG. 5: The occupations P(l) of the electronic levels I = 1, . . . 2j + 1 (see Sect. IV.) close the the 
critical point a = (3 (f?i = 0,% = 1; j • = 7/2) of the phase transition between the quasi-normal and 
radiation phase for the ground n = 1 and fourth n = 4 excited state. For f2i ^ 0,2 the occupations 
remain qualitatively the same. 
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FIG. 6: Entropies of occupation of electronic levels (|52p as function of the number of excited state 
m for j = 3/2 (4 levels). The states with entropy lower than the limiting value log(2j + 1) have 
larger measure of localization. 
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FIG. 7: Examples of the spectral density of states F{E) (|53j) for a = 2 and /? = 1.5,2.5 (We 
plot the resonant case £l\ = O2 = 1; j = 7/2; the non-resonant cases do not change the picture 
qualitatively. The same applies for previous two Figs.) 
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